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■ Single- and double-differential inclusive dijet cross sections in neutral current 

deep inelastic ep scattering have been measured with the ZEUS detector using 
an integrated luminosity of 374 pb _1 . The measurement was performed at large 
values of the photon virtuality, Q 2 , between 125 and 20 000 GeV 2 . The jets 
were reconstructed with the cluster algorithm in the Breit reference frame 
^ ■ and selected by requiring their transverse energies in the Breit frame, E^-g, to 

be larger than 8 GeV. In addition, the invariant mass of the dijet system, Mjj, 
was required to be greater than 20 GeV. The cross sections are described by the 
predictions of next-to-leading-order QCD. 



The ZEUS Collaboration 



H. Abramowicz 44 ' ^ , I. Abt 34 , L. Adamczyk 13 , M. Adamus 53 , R. Aggarwal 7 , S. Antonelli 4 , 
P. Antonioli 3 , A. Antonov 32 , M. Arneodo 49 , V. Aushev 26,aa , Y. Aushev 26 ' aa , O. Bachynska 15 , 

A. Bamberger 19 , A.N. Barakbaev 25 , G. Barbagli 17 , G. Bari 3 , F. Barreiro 29 , D. Bartsch 5 , 
M. Basile 4 , O.Behnke 15 , J. Behr 15 , U. Behrens 15 , L. Bellagamba 3 , A. Bertolin 38 , S. Bhadra 56 , 
M. Bindi 4 , C. Blohm 15 , V. Bokhonov 26 , T. Bold 13 , E.G. Boos 25 , K. Borras 15 , D. Boscherini 3 , 

D. Bot 15 , S.K. Boutle 51 , I. Brock 5 , E. Brownson 55 , R. Brugnera 39 , N. Brummer 36 , A. Bruni 3 , 

G. Bruni 3 , B. Brzozowska 52 , P.J. Bussey 20 , J.M. Butterworth 51 , B. Bylsma 36 , A.Caldwell 34 , 
M.Capua 8 , R. Carlin 39 , CD. Catterall 56 , S. Chekanov 1 , J. Chwastowski 12 ^, J. Ciborowski 52 '^', 
R. Ciesielski 15 ^, L. Cifarelli 4 , F. Cindolo 3 , A. Contin 4 , A.M. Cooper-Sarkar 37 , N. Coppola 15 '*, 
M. Corradi 3 , F. Corriveau 30 , M. Costa 48 , G. D'Agostini 42 , F. Dal Corso 38 , J. del Peso 29 , 
R.K. Dementiev 33 , S. De Pasquale 4 ^, M. Derrick 1 , R.C.E. Devenish 37 , D. Dobur 19 '", 

B. A. Dolgoshein 32 , G. Dolinska 26 , A.T. Doyle 20 , V. Drugakov 16 , L.S. Durkin 36 , S. Dusini 38 , 
Y. Eisenberg 54 , P.F. Ermolov 33 ^, A. Eskreys 12 , S. Fang 15j , S. Fazio 8 , J. Ferrando 37 , 
M.I. Ferrero 48 , J. Figiel 12 , M.Forrest 20 , B.Foster 37 , S. Fourletov 50 '™, G. Gach 13 , A. Galas 12 , 

E. Gallo 17 , A. Garfagnini 39 , A. Geiser 15 , I. Gialas 21 ' x , L.K. Gladilin 33 , D. Gladkov 32 , 

C. Glasman 29 , O. Gogota 26 , Yu.A. Golubkov 33 , P. G6ttlicher 15 ' fc , I. Grabowska-Bold 13 , 
J. Grebenyuk 15 , I. Gregor 15 , G. Grigorescu 35 , G. Grzelak 52 , O. Gueta 44 , C. Gwenlan 37 ' ac , 
T. Haas 15 , W. Hain 15 , R. Hamatsu 47 , J.C. Hart 43 , H. Hartmann 5 , G. Hartner 56 , E. Hilger 5 , 

D. Hochman 54 , R. Hori 46 , K. Horton 37 ' ad , A. Hiittmann 15 , G. Iacobucci 3 , Z.A. Ibrahim 10 , 
Y. Iga 41 , R. Ingbir 44 , M. Ishitsuka 45 , H. -P. Jakob 5 , F. Januschek 15 , M.Jimenez 29 , T.W.Jones 51 , 
M. Jiingst 5 , I. Kadenko 26 , B. Kahle 15 , B. Kamaluddin 10 't, S. Kananov 44 , T. Kanno 45 , 
U. Karshon 54 , F. Karstens 19 ^, I.I. Katkov 15 -', M. Kaur 7 , P. Kaur 7,d , A. Keramidas 35 , 
L.A. Khein 33 , JY. Kim 9 , D. Kisielewska 13 , S. Kitamura 47 ' a/ \ R. Klanner 22 , U. Klein 15 -" 1 , 

E. Koffeman 35 , P. Kooijman 35 , Ie. Korol 26 , LA. Korzhavina 33 , A. Kotanski 14,9 , U. Kotz 15 , 

H. Kowalski 15 , P. Kulinski 52 , O. Kuprash 26 ' afc , M. Kuze 45 , A. Lee 36 , B.B. Levchenko 33 , 
A. Levy 44 , V. Libov 15 , S. Limentani 39 , TY. Ling 36 , M. Lisovyi 15 , E. Lobodzinska 15 , 
W. Lohmann 16 , B. Lohr 15 , E. Lohrmann 22 , J.H. Loizides 51 , K.R. Long 23 , A. Longhin 38 , 

D. Lontkovskyi 26 > a \ OYu. Lukina 33 , P. Luzniak 52 ' afe , J. Maeda 45 ' 0-9 , S. Magill 1 , I. Makarenko 26,ab , 

J. Malka 52 ' afc , R. Mankel 15 , A. Margotti 3 , G. Marini 42 , J.F. Martin 50 , A. Mastroberardino 8 , 

M.C.K. Mattingly 2 , I.-A. Melzer-Pellmann 15 , S. Miglioranzi 15 '", F. Mohamad Idris 10 , V.Monaco 48 , 

A. Montanari 15 , J.D. Morris 6 ' c , K. Mujkic 15 ' , B. Musgrave 1 , K. Nagano 24 , T. Namsoo 15 ' p , 

R. Nania 3 , D. Nicholass 1 ' , A. Nigro 42 , Y. Ning 11 , U. Noor 56 , D. Notz 15 , R.J. Nowak 52 , 

A.E. Nuncio-Quiroz 5 , BY. Oh 40 , N. Okazaki 46 , K.Oliver 37 , K. Olkiewicz 12 , Yu. Onishchuk 26 , 

K. Papageorgiu 21 , A. Parenti 15 , E. Paul 5 , J.M. Pawlak 52 , B. Pawlik 12 , P. G. Pelfer 18 , 

A. Pellegrino 35 , W. Perlanski 52 ' afc , H. Perrey 22 , K. Piotrzkowski 28 , P. Plucinski 53 ' a ', N.S. Pokrovskiy 25 , 

A. Polini 3 , A.S. Proskuryakov 33 , M. Przybycieh 13 , A. Raval 15 , D.D. Reeder 55 , B. Reisert 34 , 

Z. Ren 11 , J. Repond 1 , Y.D. Ri 47 ' ai , A. Robertson 37 , P. Roloff 15 , E. Ron 29 , I. Rubinsky 15 , 

M. Ruspa 49 , R. Sacchi 48 , A. Salii 26 , U. Samson 5 , G. Sartorelli 4 , A. A. Savin 55 , D.H. Saxon 20 , 

M. Schioppa 8 , S. Schlenstedt 16 , P. Schleper 22 , W.B. Schmidke 34 , U. Schneekloth 15 , V. Schonberg 5 , 

T. Schorner-Sadenius 15 , J. Schwartz 30 , F. Sciulli 11 , L.M. Shcheglova 33 , R. Shehzadi 5 , 



I 



S. Shimizu 46 '™, I. Singh 7 ' d , 1.0. Skillicorn 20 , W. Slominski 14 , W.H. Smith 55 , V. Sola 48 , 
A. Solano 48 , D. Son 27 , V. Sosnovtsev 32 , A. Spiridonov 15 ' 9 , H. Stadie 22 , L. Stanco 38 , 
A. Stern 44 , T.P.Stewart 50 , A. Stifutkin 32 , P. Stopa 12 , S.Suchkov 32 , G. Susinno 8 , L. Suszycki 13 , 
J. Sztuk-Dambietz 22 , D. Szuba 15 ' r , J. Szuba 15 - 5 , A.D. Tapper 23 , E. Tassi 8 ' 6 , J. Terron 29 , 
T. Theedt 15 , H. Tiecke 35 , K. Tokushuku 24 ^, O.Tomalak 26 , J. Tomaszewska 15 '*, T. Tsurugai 31 , 
M. Turcato 22 , T. Tymieniecka 53 ' am , C. Uribe-Estrada 29 , M. Vazquez 35 '", A. Verbytskyi 15 , 
O. Viazlo 26 , N.N. Vlasov 19 ^, O. Volynets 26 , R. Walczak 37 , W.A.T. Wan Abdullah 10 , 
J.J. Whitmore 40 ' ae , J. Whyte 56 , L. Wiggers 35 , M. Wing 51 , M. Wlasenko 5 , G. Wolf 15 , 
H.Wolfe 55 , K. Wrona 15 , A.G. Yagiies-Molina 15 , S. Yamada 24 , Y. Yamazaki 24 ' 2 , R. Yoshida 1 , 
C. Youngman 15 , A.F. Zarnecki 52 , L. Zawiejski 12 , O. Zenaiev 26 , W. Zeuner 15 '™, B.O. Zhautykov 25 , 
N. Zhmak 26 < aa , C. Zhou 30 , A. Zichichi 4 , M. Zolko 26 , D.S. Zotkin 33 , Z. Zulkapli 10 



II 



Argonne National Laboratory, Argonne, Illinois 60439-4815, USA 

Andrews University, Berrien Springs, Michigan 49104-0380, USA 

INFN Bologna, Bologna, Italy B 

University and INFN Bologna, Bologna, Italy B 

Physikalisches Institut der Universitdt Bonn, Bonn, Germany c 

H.H. Wills Physics Laboratory, University of Bristol, Bristol, United Kingdom D 

Panjab University, Department of Physics, Chandigarh, India 

Calabria University, Physics Department and INFN, Cosenza, Italy B 

Institute for Universe and Elementary Particles, Chonnam National University, 

Kwangju, South Korea 

Jabatan Fizik, Universiti Malaya, 50603 Kuala Lumpur, Malaysia E 

Nevis Laboratories, Columbia University, Irvington on Hudson, New York 10027, 

USA F 

The Henryk Niewodniczanski Institute of Nuclear Physics, Polish Academy of 
Sciences, Cracow, Poland G 

Faculty of Physics and Applied Computer Science, AGH-University of Science and 
Technology, Cracow, Poland H 

Department of Physics, Jagellonian University, Cracow, Poland 
Deutsches Elektronen- Synchrotron DESY, Hamburg, Germany 
Deutsches Elektronen- Synchrotron DESY, Zeuthen, Germany 
INFN Florence, Florence, Italy B 
University and INFN Florence, Florence, Italy B 

Fakultat fur Physik der Universitdt Freiburg i.Br., Freiburg i.Br., Germany 
School of Physics and Astronomy, University of Glasgow, Glasgow, United King- 
dom D 

Department of Engineering in Management and Finance, Univ. of the Aegean, 
Chios, Greece 

Hamburg University, Institute of Experimental Physics, Hamburg, Germany 1 
Imperial College London, High Energy Nuclear Physics Group, London, United 
Kingdom D 

Institute of Particle and Nuclear Studies, KEK, Tsukuba, Japan J 
Institute of Physics and Technology of Ministry of Education and Science of Kaza- 
khstan, Almaty, Kazakhstan 

Institute for Nuclear Research, National Academy of Sciences, and Kiev National 
University, Kiev, Ukraine 

Kyungpook National University, Center for High Energy Physics, Daegu, South Ko- 

K 

rea 

Institut de Physique Nucleaire, Universite Catholique de Louvain, Louvain-la-Neuve, 
Belgium L 

Departamento de Fisica Teorica, Universidad Autonoma de Madrid, Madrid, 
Spain M 

Department of Physics, McGill University, Montreal, Quebec, Canada H3A 2T8 N 
Meiji Gakuin University, Faculty of General Education, Yokohama, Japan J 
Moscow Engineering Physics Institute, Moscow, Russia 



III 



Moscow State University, Institute of Nuclear Physics, Moscow, Russia p 
Max- Planck- Institut fur Physik, Miinchen, Germany 
NIKHEF and University of Amsterdam, Amsterdam, Netherlands ® 
Physics Department, Ohio State University, Columbus, Ohio 43210, USA A 
Department of Physics, University of Oxford, Oxford, United Kingdom 
INFN Padova, Padova, Italy B 

Dipartimento di Fisica dell' Universita and INFN, Padova, Italy B 

Department of Physics, Pennsylvania State University, University Park, 

Pennsylvania 16802, USA F 

Polytechnic University, Sagamihara, Japan J 

Dipartimento di Fisica, Universita 'La Sapienza' and INFN, Rome, Italy 
Rutherford Appleton Laboratory, Chilton, Didcot, Oxon, United Kingdom D 
Raymond and Beverly Sackler Faculty of Exact Sciences, School of Physics, 
Tel Aviv University, Tel Aviv, Israel R 

Department of Physics, Tokyo Institute of Technology, Tokyo, Japan J 
Department of Physics, University of Tokyo, Tokyo, Japan J 
Tokyo Metropolitan University, Department of Physics, Tokyo, Japan J 
Universita di Torino and INFN, Torino, Italy B 

Universita del Piemonte Orientale, Novara, and INFN, Torino, Italy B 
Department of Physics, University of Toronto, Toronto, Ontario, Canada M5S 
1A7 N 

Physics and Astronomy Department, University College London, London, United 
Kingdom D 

Warsaw University, Institute of Experimental Physics, Warsaw, Poland 
Institute for Nuclear Studies, Warsaw, Poland 

Department of Particle Physics, Weizmann Institute, Rehovot, Israel s 
Department of Physics, University of Wisconsin, Madison, Wisconsin 53706, USA A 
Department of Physics, York University, Ontario, Canada M3J IPS N 



IV 



supported by the US Department of Energy 

B supported by the Italian National Institute for Nuclear Physics (INFN) 

c supported by the German Federal Ministry for Education and Research (BMBF), 

under contract No. 05 H09PDF 
D supported by the Science and Technology Facilities Council, UK 

E supported by an FRGS grant from the Malaysian government 

F supported by the US National Science Foundation. Any opinion, findings and con- 
clusions or recommendations expressed in this material are those of the authors and 

do not necessarily reflect the views of the National Science Foundation. 
G supported by the Polish Ministry of Science and Higher Education as a scientific 

project No. DPN/N188/DESY/2009 
H supported by the Polish Ministry of Science and Higher Education as a scientific 
project (2009-2010) 

1 supported by the German Federal Ministry for Education and Research (BMBF), 
under contract No. 05h09GUF, and the SFB 676 of the Deutsche Forschungsge- 
meinschaft (DFG) 

J supported by the Japanese Ministry of Education, Culture, Sports, Science and 

Technology (MEXT) and its grants for Scientific Research 
K supported by the Korean Ministry of Education and Korea Science and Engineering 

Foundation 

supported by FNRS and its associated funds (IISN and FRIA) and by an Inter- 
University Attraction Poles Programme subsidised by the Belgian Federal Science 
Policy Office 

M supported by the Spanish Ministry of Education and Science through funds provided 
by CICYT 

N supported by the Natural Sciences and Engineering Research Council of Canada 
(NSERC) 

partially supported by the German Federal Ministry for Education and Research 
(BMBF) 

p supported by RF Presidential grant N 41-42.2010.2 for the Leading Scientific Schools 

and by the Russian Ministry of Education and Science through its grant for Scientific 

Research on High Energy Physics 
® supported by the Netherlands Foundation for Research on Matter (FOM) 

R supported by the Israel Science Foundation 

s supported in part by the MINERVA Gesellschaft fur Forschung GmbH, the Israel 
Science Foundation (grant No. 293/02-11.2) and the US-Israel Binational Science 
Foundation 



V 



a also affiliated with University College London, United Kingdom 

h now at University of Salerno, Italy 

c now at Queen Mary University of London, United Kingdom 

d also working at Max Planck Institute, Munich, Germany 

e also Senior Alexander von Humboldt Research Fellow at Hamburg University, Insti- 
tute of Experimental Physics, Hamburg, Germany 

* also at Cracow University of Technology, Faculty of Physics, Mathemathics and 

Applied Computer Science, Poland 

9 supported by the research grant No. 1 P03B 04529 (2005-2008) 

h now at Rockefeller University, New York, NY 10065, USA 

i now at DESY group FS-CFEL-1 

J now at Institute of High Energy Physics, Beijing, China 

k now at DESY group FEB, Hamburg, Germany 

1 also at Moscow State University, Russia 

m now at University of Liverpool, United Kingdom 

" now at CERN, Geneva, Switzerland 

° also affiliated with Universtiy College London, UK 

p now at Goldman Sachs, London, UK 

q also at Institute of Theoretical and Experimental Physics, Moscow, Russia 

r also at INP, Cracow, Poland 

s also at FPACS, AGH-UST, Cracow, Poland 

' partially supported by Warsaw University, Poland 

" now at Istituto Nucleare di Fisica Nazionale (INFN), Pisa, Italy 

v now at Haase Energie Technik AG, Neumiinster, Germany 

w now at Department of Physics, University of Bonn, Germany 

x also affiliated with DESY, Germany 

y also at University of Tokyo, Japan 

z now at Kobe University, Japan 

' deceased 

aa supported by DESY, Germany 

ab supported by the Bogolyubov Institute for Theoretical Physics of the National 

Academy of Sciences, Ukraine 

ac STFC Advanced Fellow 

ad nee Korcsak-Gorzo 

ae This material was based on work supported by the National Science Foundation, 

while working at the Foundation. 

a f also at Max Planck Institute, Munich, Germany, Alexander von Humboldt Research 

Award 

ag now at Tokyo Metropolitan University, Japan 

ah now at Nihon Institute of Medical Science, Japan 

a% now at Osaka University, Osaka, Japan 

aj also at Lodz University, Poland 



VI 



* member of Lodz University, Poland 

1 now at Lund University, Lund, Sweden 

n also at University of Podlasie, Siedlce, Poland 



VII 



1 Introduction 



Measurements of jet cross sections are a well established tool for QCD studies and have 
been performed for many different observables at HERA [THTO]. For jet analyses in neutral 
current (NC) deep inelastic scattering (DIS), the Breit reference frame [IT] is preferred, 
since it provides a maximal separation between the products of the beam fragmentation 
and the hard jets. In this frame, the exchanged space-like virtual boson has 3-momentum 
q = (0, 0, —Q) and is collinear with the incoming parton. While retaining hard QCD 
processes at leading order in the strong coupling constant a s , the contribution from the 
parton-model process can be suppressed by requiring the production of jets with high 
transverse energy in the Breit frame. Therefore, measurements of jet cross sections in the 
Breit frame are directly sensitive to hard QCD processes, allowing tests of perturbative 
QCD (pQCD), of the factorisation ansatz and of the parton distribution functions (PDFs) 
of the proton. 

At large boson virtualities, Q 2 , the experimental and theoretical systematic uncertainties 
are small and, thus, use of the large HERA data sample can provide powerful physical 
constraints. Jet cross-section data from the high-Q 2 region have been included in the 
ZEUS PDF fit, thereby significantly reducing the uncertainty on the gluon density in the 
medium- to high-x region [1'2\ . 

Measurements of dijet production in DIS at HERA have so far been performed with either 
smaller data sets [HEj or normalised to the inclusive NC DIS cross section [SJ. In this 
paper, measurements of inclusive dijet production at large values of Q 2 are presented 
using an integrated luminosity of 374 pb _1 . Here, differential dijet cross sections as a 
function of Q 2 , of the mean jet transverse energy of the dijet system in the Breit frame, 
-Ey B , of the dijet invariant mass, My, of the half- difference of the jet pseudorapidities in 

the Breit frame, r/* = tf^ 1 — -rf^ 2 /2, of the fraction of the proton momentum taken by 

the interacting parton, £ = xej (l + (Mjj) 2 /<5 2 ), and of xbj are presented. Here, xbj is 
the Bjorken scaling variable that defines, for the parton-model process, the fraction of 
the proton momentum carried by the struck massless parton. Measurements of the dijet 
cross section as a function of £ and E^ B are also shown for different regions of Q 2 . 

The measured single- and double-differential cross sections are compared with next-to- 
leading-order (NLO) QCD calculations. 



2 Experimental set-up 

A detailed description of the ZEUS detector can be found elsewhere [13J. A brief outline 
of the components that are most relevant for this analysis is given below. 
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Charged particles were tracked in the central tracking detector (CTD) [H], which operated 
in a magnetic field of 1.43 T provided by a thin superconducting solenoid. The CTD 
consisted of 72 cylindrical drift-chamber layers, organised in nine superlayers covering the 
polar-angle region 15° < < 164°. For data taken during the years 1998 to 2000, tracks 
were reconstructed using the CTD only. Starting from the year 2004, the CTD and a 
silicon microvertex detector (MVD) [15], installed between the beam-pipe and the inner 
radius of the CTD, were used. 

The high-resolution uranium-scintillator calorimeter (CAL) [IB] covered 99.7% of the 
total solid angle and consisted of three parts: the forward (FCAL), the barrel (BCAL) 
and the rear (RCAL) calorimeters. Each part was subdivided transversely into towers 
and longitudinally into one electromagnetic section (EMC) and either one (in RCAL) 
or two (in BCAL and FCAL) hadronic sections (HAC). The smallest subdivision of the 
calorimeter is called a cell. Under test-beam conditions, the CAL single-particle relative 
energy resolutions were a(E)/E = 0.18/\^E for (anti-)electrons and a(E)/E = 0.35/\/E 
for hadrons, with E in GeV. 

The luminosity was measured using the Bethe- Heftier reaction ep — > e'yp by the luminosity 
detector which consisted of a lead-scintillator [17] calorimeter and, in the 2004-2007 
running period, an independent magnetic spectrometer [18]. 

The electronic beam in HERA was naturally transversely polarised through the Sokolov- 
Ternov effect [T9~||2"0"]. The characteristic build-up time for the HERA accelerator was 
approximately 40 minutes. Starting from the year 2004, spin rotators on either side of 
the ZEUS detector changed the transverse polarisation of the beam into longitudinal 
polarisation and back again. The electron beam polarisation was measured using two 
independent polarimeters, the transverse polarimeter (TPOL) [21] and the longitudinal 
polarimeter (LPOL) [22] • Both devices exploited the spin-dependent cross section for 
Compton scattering of circularly polarised photons off electrons to measure the beam po- 
larisation. The luminosity and polarisation measurements were made over time intervals 
that were much shorter than the polarisation build-up time. 



3 Event selection and reconstruction 

The data used in this analysis were collected during the periods 1998-2000 and 2004-2007, 
when HERA operated with protons of energy E p = 920 GeV and electrons or positrons 
of energy E e = 27.5 GeV, and correspond to an integrated luminosity of 203 pb" 1 for the 

1 Here and in the following, the term "electron" denotes generically both the electron and the positron 
unless otherwise stated. 
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electron data and 171 pb 1 for the positron data. The mean luminosity-weighted average 
polarisation of the data was —0.03. 

A three- level trigger system was used to select events online [13l[23]. At the third level, 
NC DIS events were accepted on the basis of the identification of a scattered-electron 
candidate using localised energy depositions in the CAL. At the second level, charged- 
particle tracks were reconstructed online by using the ZEUS global tracking trigger [21], 
which combined information from the CTD and MVD. These online tracks were used to 
reconstruct the interaction vertex and reject non-ep background. At the first level, only 
coarse calorimeter and tracking information was available. Events were selected using 
criteria based on the energy and transverse energy measured in the CAL. Starting from 
the year 2004, additional tracking requirements were introduced to adapt the trigger rates 
to the higher instantaneous luminosity. 

Events were selected offline using criteria which were slightly changed with respect to 
those used in the previous ZEUS dijet measurement [6]. These changes reflect the new 
phase-space definition of the measurements adopted here. The main steps of the selection 
are briefly listed below. 

The scattered-electron candidate was identified from the pattern of energy deposits in the 
CAL [25] . The energy, E' e , thepolaii angle, 9 e , and the azimuthal angle, e , of the electron 
candidate were determined from the CAL measurements. For events in which the electron 
was found inside the CTD acceptance, the angles 9 e and e were reconstructed from the 
associated electron track. The photon virtuality Q 2 and the Bjorken scaling variable XBj 
were reconstructed using the double angle (DA) method [26J. The inelasticity variable y 
was determined from the condition y = Q 2 /xb]S, where s is the square of the centre-of- 
mass energy. 

3.1 Inclusive event selection 

The phase space of the measurement was 125 < Q 2 < 20 000 GeV 2 and 0.2 < y < 0.6. 
Events were selected if: 

• an electron candidate of energy E' e > 10 GeV was found. This requirement ensured 
a high and well understood electron-finding efficiency and suppressed background 
from photoproduction events in which the scattered electron escapes down the rear 
beampipe; 

2 The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton-beam direction, referred to as the "forward direction" , and the X axis pointing left towards 
the center of HERA. The coordinate origin is at the nominal interaction point. 
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• the total energy not associated with the electron candidate within a cone of radius 0.7 
units in the pseudorapidity-azimuth (rj — tp) plane around the electron direction was 
less than 10% of the electron energy. This condition removed photoproduction and 
DIS events in which part of a jet was falsely identified as the scattered electron; 

• a track matched to the energy deposit in the CAL was found in events in which an 
electron was reconstructed within the acceptance of the tracking detectors. This was 
done by restricting the distance of closest approach between the track extrapolated 
to the CAL surface and the energy cluster position to within 10 cm, and requiring an 
electron track momentum greater than 3 GeV; 

• the vertex position along the beam axis was in the range that was given by the nom- 
inal vertex position plus or minus three times the width of the vertex distribution 
approximated by a Gaussian. Both the nominal vertex position and the width of 
the distribution varied between the different data-taking periods. Typical values were 
\Zvtx\ < 30 cm. This condition helped to select events consistent with ep interactions; 

• -Pr.miss/v^T < 2.5 GeV 1 / 2 , where -Pr.miss is the missing transverse momentum as 
measured with the CAL and Et is the total transverse energy in the CAL. This cut 
removed cosmic-ray events and beam-related background; 

• 38 < (E — Pz) < 65 GeV, where E is the total energy, E = "Y^i^ii an ^ Pz is the 
Z-component of the vector P = J^Pi- The sums run over all clusters of energy 
deposits in the CAL. This cut removed events with large initial-state radiation and 
further reduced the background from photoproduction. 

In addition, events were rejected if a second electron candidate with azimuthal separation 
A0 > 3 from the first candidate was found, a ratio of transverse momenta of the two 
candidates between 0.8 and 1.2 was measured, and, in addition, the rest of the CAL 
energy, besides the two electromagnetic energy clusters, was below 3 GeV. This condition 
removed elastic Compton scattering events (ep — > ejp). 

3.2 Jet selection 

The kj* cluster algorithm [27] was used in the longitudinally invariant inclusive mode [28] 
to reconstruct jets in the hadronic final state assuming massless objects. In data, the 
algorithm was applied to the CAL cells after excluding those associated with the scattered- 
electron candidate. The jet search was performed in the r] — ip plane of the Breit frame. 
The jet variables were defined according to the Snowmass convention [29] . 

After reconstructing the jet variables in the Breit frame, the massless four-momenta were 
boosted into the laboratory frame, where the transverse energy, E^ LAB , and the pseudo- 
rapidity, ?7l2b) °f eacri J e t were calculated. Energy corrections [2"}I3T)|I3"T] were then applied 



4 



to the jets in the laboratory frame and propagated into E^ B , the transverse jet energy in 
the Breit frame, in order to compensate for energy losses in the inactive material in front 
of the CAL. 

The following cuts were applied to select a clean sample of high-Q 2 DIS jet events: 

• events were removed from the sample if the distance AR of any of the jets to the 
electron candidate in the r\ — tp plane of the laboratory frame was smaller than 1 unit, 

AR = J (t/lab — V e ) 2 + (fhAB ~ P e ) 2 < 1> where (p e and rf are the azimuthal angle 
and pseudorapidity of the scattered-electron candidate, respectively This require- 
ment removed some background from photoproduction and improved the purity of 
the sample; 

• events were removed from the sample if a jet was in the backward region of the detector 
(Vlab < ~~ -0- This requirement removed events in which a radiated photon from the 
electron was misidentified as a jet in the Breit frame; 

• Ej, lab > ^ GeV. This cut removed a small number of jets for which the uncertainty 
on the energy correction was large. 

The dijet sample was then selected requiring the events to fullfil the following conditions 
(which also define the dijet phase space of the measurement): 

• at least two jets in the pseudorapidity range —1 < ^lab < 2-5 were found; 

• of these jets, the two with the highest transverse energies in the Breit frame were 
required to have Ej? B > 8 GeV; 

• the invariant dijet mass, Mjj, of the two highest-transverse-energy jets in the event was 
required to exceed 20 GeV. This requirement was introduced to make the theoretical 
calculations infrared insensitive. Despite this cut the NLO QCD calculations in the 
region 0.65 < 77* < 2 still exhibited infrared sensitivity (Section 

The final sample consisted of 31440 dijet events. 

4 Monte Carlo simulations and data corrections 

Samples of Monte Carlo (MC) events were generated to determine the response of the 
detector to jets of hadrons and the correction factors necessary to obtain the hadron-level 
jet cross sections. The hadron level is defined in terms of hadrons with lifetime r > 10 ps. 
The generated events were passed through the GEANT-based [32] ZEUS detector- and 
trigger-simulation programs [13]. They were reconstructed and analysed by the same 
program chain as the data. 
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Neutral current DIS events including radiative effects were simulated using the Her- 
acles 4.6.1 [33] program with the Djangoh 1.1 [31] interface to the hadronisation 
programs. HERACLES includes corrections for initial- and final-state radiation, vertex 
and propagator terms, and two-boson exchange. The QCD cascade is simulated using the 
colour-dipole model (CDM) [32] including the LO QCD diagrams as implemented in Ari- 
adne 4.08 [36] and, alternatively, with the matrix-element plus parton-shower (MEPS) 
approach of Lepto 6.5 [31] • The CTEQ5D [3H] proton PDFs were used for these sim- 
ulations. Fragmentation into hadrons is performed using the Lund string model [39] as 
implemented in Jetset 7.4 [40] . 

The jet search was performed on the MC events using the energy measured in the CAL 
cells in the same way as for the data. The same jet algorithm was also applied to the 
final-state particles (hadron level) and to the partons available after the parton shower 
simulation (parton level). 

The data were corrected to the hadron level and for QED-radiative effects using bin-by- 
bin correction factors obtained from the MC samples. For this approach to be valid, the 
uncorrected distributions of the data must be well described by the MC simulations. This 
condition was in general satisfied by both the Ariadne and Lepto MC. Figures [T]|3] 
show comparisons of data with MC simulations for all observables in which cross sections 
are presented in this paper. The MC simulations give a good description of all the data 
distributions. The Lepto model gives a slightly better description of the data and was 
thus used as the default model; Ariadne was then used to estimate the systematic effect 
on the correction procedure due to the parton-shower model. In all cases, the correction 
factors differed from unity by (5 — 30) %. These correction factors took into account the 
efficiency of the trigger, the selection criteria and the purity and efficiency of the jet 
reconstruction. The QED corrections typically amounted to between 3 and 6%. 

Starting in 2004, HERA provided longitudinally polarised lepton beams. The effect of the 
polarisation on the measured data events was corrected for by adjusting the data event 
weights with the ratio of the predictions for the unpolarised and polarised cross sections 
as determined with the Hector program |41j . 

5 NLO QCD calculations 

Next-to-leading-order (0(a 2 )) QCD calculations were obtained using the program Nlo- 
jetH — h [12] . The calculations were performed in the MS renormalisation and factorisation 
schemes. The number of flavours was set to five and the factorisation scale was chosen 
to be fip = Q. Calculations with different choices of the renormalisation scale, /j,r, were 

performed: the default choice was ji 2 R = Q 2 + E^ B . Alternatively, the scales Q 2 and 
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E^ B were investigated. The strong coupling constant was calculated at two loops with 
A|jgr = 226 MeV, corresponding to a s (M z ) = 0.118. The calculations were performed 
using the CTEQ6.6 [13] parameterisations of the proton PDFs. The kx cluster algorithm 
was also applied to the partons in the events generated by NlojetH — h in order to obtain 
the jet cross-section predictions. The predictions of NlojetH — h were cross-checked with 
the DlSENT program [33]. Both programs agreed to better than 1%. 

The lack of sensitivity to the infared cutoff of the NLO QCD calculations was verified 
by determining the total inclusive dijet cross section as a function of the My cut in each 
analysis bin separately. Except for the two highest rj* bins, the theoretical predictions 
were found to be infrared insensitive [15]. 

The data presented in this paper are, among others, intended for the use in QCD PDF 
fits, aiming specifically at a further improvement of the uncertainty on the gluon density 
at large values of x. In order to demonstrate the usefulness of the data for this purpose, 
Figs. HI and [5] show, as a function of the variables £ and -Ej^ B in different regions of 
Q 2 , the NLO predictions for the fraction of events which are initiated by a gluon from 
the proton using the CTEQ6.6 PDFs. This gluon fraction ranges from about 75% at 
125 < Q 2 < 250 GeV 2 and small £ to about 5% at the highest Q 2 above 5000 GeV 2 , where 
£ is approximately confined to values above 0.1. In the lower Q 2 regions, the gluon fraction 
is also significant for large values of £. Since this region is not statistically limited, precise 
input for the PDF fits can be expected. Figure [6] shows the relative CTEQ6.6 PDF 
uncertainty, the uncertainty due to missing higher orders in the calculation estimated by 
variation of fi R and the theoretical predictions from MSTW2008 [3H], ZEUS- JETS [TJ 
and ZEUS-S [H]. The corresponding uncertainties for the latter three PDF sets are not 
shown. The CTEQ6.6 PDF uncertainty and the observed spread between the various 
PDF sets is in some regions of the considered dijet phase space larger than the uncertainty 
arising from missing higher orders. 

The measurements refer to jets of hadrons, whereas the NLO QCD calculations refer to 
jets of partons. The QCD predictions were corrected to the hadron level using the MC 
hadronisation model described in the previous section to give multiplicative factors, Ch a dr, 
defined as the ratio of the cross section for jets of hadrons over that for jets of partons. 
The ratios obtained with Ariadne and Lepto were averaged to obtain the Ch a d r factors, 
which differ from unity by less than 5 %. 

Neither NlojetH — \- nor Disent includes the contributions from Z° exchange; MC simu- 
lated events with and without Z° exchange were used to include this effect in the pQCD 
predictions. In the following, pQCD calculations refer to the fully corrected predictions, 
unless otherwise stated. 

Several sources of uncertainty in the theoretical predictions were considered: 
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• the uncertainty on the NLO QCD calculations due to terms beyond NLO, estimated 
by varying [Ir by a factor of two up and down, was below ±6% at low Q 2 and low 
Ej, B and decreased to below ±3% in the highest-Q 2 region; 

• the uncertainty on the NLO QCD calculations due to that on a s (Mz) was estimated 
by repeating the calculations using two additional sets of proton PDFs, CTEQ6.6A3 
and CTEQ6.6A2, determined assuming a s (Mz) = 0.114 and 0.122, respectively. The 
difference between the calculations using these sets and CTEQ6.6 was scaled to reflect 
the current uncertainty on a s jlH]. The resulting uncertainty on the cross sections was 
mostly below ±3%; 

• the uncertainty on the modelling of the parton shower was estimated as half the 
difference between the correction factors calculated from the Lepto and Ariadne 
models. The resulting uncertainty on the cross sections was typically less than 2%; 

• the uncertainty on the NLO calculations due to the proton PDFs was estimated by 
repeating the calculations using 44 additional sets from the CTEQ6.6 analysis, which 
takes into account the statistical and correlated systematic experimental uncertainties 
of each data set used in the determination of the proton PDFs. The resulting uncer- 
tainty on the cross sections was about ±4% at low Q 2 and decreased to around ±2% 
at high Q 2 ; 

• the uncertainty of the calculations in the value of \xp was estimated by repeating the 
calculations with fip = Q/2 and 2Q. The effect on the calculations was negligible. 

The total theoretical uncertainty was obtained by adding in quadrature the individual 
uncertainties listed above. 

6 Experimental uncertainties 

The following sources of systematic uncertainty were considered for the measured cross 
sections [45] : 

• the uncertainty on the absolute energy scale of the jets was estimated to be ±1% 
for E^ LAB > 10 GeV and ±3% for lower E^ LAB values OEMES]. The resulting 
uncertainty on the cross sections was about ±4% and increased to approximately ±6% 
in certain regions of the dijet phase space; 

• the uncertainty in the absolute energy scale of the electron candidate was estimated 
to be ±1% |SU] (±2% [55) for the data from the years 1998-2000 (2004-2007). The 
resulting uncertainty was below ±1%; 

• the differences in the results obtained by using either Ariadne or Lepto to correct 
the data for detector effects were typically below ±5%; 
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• the analysis was repeated using an alternative technique [52] to select the scattered- 
electron candidate. The resulting uncertainty was typically below ±1%; 

• the lab cut was changed to 2 and 4 GeV. The resulting uncertainty was mostly 
smaller than ±1%; 

• the uncertainty due to the selection cuts was estimated by varying the values of the 
cuts within the resolution of each variable. The effect on the cross sections was in 
general below ±2%; 

• the combined, luminosity-weighted systematic error on the polarisation measurement 
was 3.9%. The effect on the cross sections was negligible; 

• the simulation of the first-level trigger was corrected in order to match the measured 
efficiency in the data. The systematic effect on the cross sections was typically less 
than 1%. 

The systematic uncertainties not associated with the absolute energy scale of the jets were 
added in quadrature. Figure [7| shows the statistical uncertainty, the correlated systematic 
uncertainty which is caused by the jet energy scale and the quadratic sum of the correlated 
and uncorrelated systematic uncertainties as a function of Q 2 . Except for the high-Q 2 
region, the correlated uncertainty was the dominating contribution to the total experi- 
mental uncertainty. In addition, there was an overall normalisation uncertainty of ±2.2% 
for the 1998-2000 data and of ±2.6% for the 2004-2007 data. Therefore, the combined, 
luminosity-weighted average systematic uncertainty on the luminosity measurement was 
±2.5%, which was not included in the cross-section figures or the tables. 

7 Results 

The differential inclusive dijet cross sections were measured in the kinematic region 125 < 
Q 2 < 20 000 GeV 2 and 0.2 < y < 0.6. The jets were reconstructed using the cluster 
algorithm in the longitudinally invariant inclusive mode and the cross sections refer to 
jets with E^ B > 8 GeV and —1 < 77lab < 2.5. The invariant dijet mass of the two 
highest-transverse-energy jets in the event was required to be greater than 20 GeV. These 
cross sections were corrected for detector and QED radiative effects and the running of 
o em . 

7.1 Single-differential dijet cross sections 

The measurements of the single-differential inclusive dijet cross sections are presented 
in Figs El to [TU] and Tables [T] to [6] as functions of several kinematic and dijet variables. 
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Single-differential cross sections are shown for Q 2 , xsj, the mean transverse jet energy in 
the Breit frame of the two jets, E^ B , the dijet invariant mass, Mjj, the half- difference 
of the jet pseudorapidities in the Breit frame, 77*, and the logarithm of the variable £. 
The data are compared to NLO QCD calculations. The relative differences between the 
measured differential cross sections and the NLO QCD calculations are also shown. 

The single-differential dijet cross-sections da/dQ 2 and da/dx-Qj are shown in FigsEk and 
[8b. The cross-section da/dQ 2 has total experimental systematic uncertainties of the order 
of 5% (7%) at low (high) values of Q 2 . The total theoretical uncertainty is of the order 
of 7% (4%) at low (high) Q 2 . 

For the cross-section da/dx-B], most of the data points have experimental uncertainties of 
less than 5%, and also the precision of the theory predictions is better than 5% over most 
of the xb] range. 

Figures [Hk and|Hb show the single-differential dijet cross-sections da/dE 1 ^ and da/dM^. 
These measurements are particularly well suited for testing the matrix elements in the 
perturbative calculations. Mean transverse jet energies -E^b (dijet invariant masses Mjj) 
of up to 60 GeV (120 GeV) are reached with this measurement. At the largest values of 
E^ B (Mjj), experimental uncertainties of 8% (5%) are achieved; for smaller values, the 
uncertainties are even smaller. The theoretical uncertainties are approximately constant 
over the range studied and are of the order of (5 — 7) %. 

The differential dijet cross-section as a function of 77* is shown in Fig.fTOk. The experimen- 
tal uncertainties are always below 5%, the total theoretical uncertainty is also typically 
around 5%. The theoretical predictions for the last two 77* bins were removed from the 
plot due to infrared sensitivity. 

The cross-section da/dlog 10 (£) (Fig. 110b) has similar uncertainties as the distributions 
described before and shows a maximum around log 10 (£) = —1.5. At lower and higher 
values, the cross section reflects the suppression by the transverse energy requirements in 
the selection and the decreasing quark and gluon densities, respectively. 

All the measured differential cross sections are well described by NLO QCD predictions. 
7.2 Double-differential dijet cross sections 

Figures [11] to [TH show the measurements of double-differential dijet cross sections as 
functions of £^' B and log 10 (£) in different Q 2 regions (see Tables [7] and [SJ . These cross 
sections will provide valuable input for the extraction of the proton PDFs. 

The log 10 (£) distributions in different Q 2 regions in Fig. [TT] show the same behaviour 
as the integrated log 10 (£) distribution in Fig. [TUb . with a distinct maximum at values 
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that increase with increasing Q 2 . The data are very precise - even in the highest Q 2 
bin from 5 000 to 20 000 GeV 2 the experimental uncertainties are between 10 and 15% 
and originate equally from the statistical and the systematical uncertainty. At lower Q 2 
values, the experimental uncertainties become as small as (2 — 3) %. Figure [12] shows 
the level of agreement between data and predictions: the theoretical uncertainties are 
typically between 5 and 10% and, within the combined uncertainties, the data are very 
well described by the theory. 

The cross sections as functions of E^ B in different regions of Q 2 , shown in Fig. [131 
fall over 2 to 3 orders of magnitude in the range considered, with a smaller slope for 
higher Q 2 values. The statistical precision of the data is between 2% at the lowest E^ B 
and Q 2 and slightly above 10% at the highest values of these variables. The systematic 
uncertainties are mostly of the order of 3 to 5%. The theoretical uncertainties (Fig. [T4"j) 
are approximately constant in E^ B ; they are of the order of 5 to 10%, with the smaller 
values at higher Q 2 . Data and theory are in good agreement over the whole measured 
range. 

8 Summary and conclusions 

Measurements of single- and double-differential cross sections for dijet production at high- 
er 2 NC DIS were made using an integrated luminosity of 374 pb -1 . The measurements 
have very small statistical and systematic uncertainties and the description of the data by 
the predictions of NLO QCD is very good, giving a powerful and stringent justification 
of the theory. These data will provide useful precision information for the determination 
of the strong coupling constant and the extraction of the proton PDFs. 
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Q 2 bin 
(GeV 2 ) 


da/dQ 2 

(pb/GeV 2 ) (Jgtat 5 S yst S E s 


Cqed 


Chadr " Cz" 


125 ... 250 
250 ... 500 
500 ... 1000 
1000 . . . 2000 
2000 . . . 5000 
5000 . . . 20000 


0.3843 ±0.0036 i°;°olo ^aoSs 

0.1193 ±0.0015 j£°gg 

03372 ±0 000S3 +0-00065 +0.00135 

U.UO^Z ICU.UUUOO -0.00065 -0.00115 

008^ ±0 001 8 +0-00010 +0.00029 

U.UUOOO 3ZU.UUU10 -0.00010 -0.00026 

001^23 ±0 43 +0-000033 +0.000030 
u.uuiozo iru.uuuu^o -0.000033 -0.000034 

000087^ ±0 0000046 +0-0°ooo 58 +0.0000014 
u.uuuuo/o iru.uuuuu^o -0.0000057 -0.0000015 


0.97 
0.95 
0.94 
0.93 
0.93 
0.92 


0.95 

0.96 
0.96 
0.98 
1.03 
1.09 



Table 1: The measured differential cross-sections da/dQ 2 for inclusive dijet 
production. The statistical, uncorrelated systematic and jet-energy-scale (ES) un- 
certainties are shown separately. The multiplicative corrections, Cqed? which have 
been applied to the data and the corrections for hadronisation and Z° effects to be 
applied to the parton-level NLO QCD calculations, Chadr ■ Czo, are shown in the 
last two columns. 



x-Q] bin 


da/dx-Q] 
(Pb) 


"stat 


^syst 


^ES 


Cqed 


Chadr ' CgO 


0.0001 .. 


. 0.01 


6580 


±54 


+44 
-45 


+351 
-317 


0.96 


0.95 


0.01 .. 
0.02 . . 


. 0.02 
. 0.035 


2229 
711 


±31 
±14 


+44 
-44 

+19 
-20 


+98 
-94 

+27 
-23 


0.94 
0.94 


0.95 
0.96 


0.035 . . 


. 0.07 


193.8 


±4.6 


+2.8 
-2.5 


+6.1 
-5.7 


0.93 


0.99 


0.07 .. 


. 0.1 


64.4 


±2.8 


+3.8 
-3.8 


+1.0 
-1.4 


0.92 


1.03 



Table 2: Inclusive dijet cross-sections da/dx^y Other details as in the caption 
to TableUl 
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fiy B bin 
(GeV) 


da / dE^ B 

(pb/GeV) <5 st at 4 y st S ES 


Cqed 


Chadr " Cz° 


8 ... 15 
15 ... 22 
22 ... 30 
30 ... 60 


10.650 ±0.083 t° Hi t 5 4 f 5 
3.595 ±0.046 ±g:8S + -oHl 
0.848 ±0.020 t H TH 
0.0896 ±0.0031 %™ t°££ 


0.95 
0.95 
0.95 
0.95 


0.95 

0.98 
0.99 
0.99 



Table 3: Inclusive dijet cross-sections da j dE^ B . Other details as in the caption 
to TableUl 



Mjj bin 
(GeV) 


da/dMjj 

(pb/GeV) 4tat 4yst $ES 


Cqed 


Chadr " CzO 


20 ... 30 
30 ... 45 
45 ... 65 
65 ... 120 


5.048 ±0.049 ±£°™ Tit 
2.693 ±0.028 t° Z t°f 2 ? 
0.726 ±0.012 t° Z 

0.0681 ±0.0020 


0.95 
0.95 
0.95 
0.95 


0.95 

0.97 
0.98 
0.97 



Table 4: Inclusive dijet cross-sections da/dMy y Other details as in the caption 
to Tabled 



rj* bin 


dajdrf 

(Pb) <5 st at 4yst 5es 


Cqed 


Chadr " CzO 


... 0.2 
0.2 ... 0.4 
0.4 ... 0.65 
0.65 ... 0.95 
0.95 ... 2 


106.1 ±1.6 ^ tH 
105.4 ±1.6 It? 
101.0 ±1.4 tto 
78.2 ±1.1 «i ±H 
17 14 ±0 2 7 +0 - 42 + 1 - 06 

Lt.L<± ZCV.ZI „ 42 _ 1Q2 


0.95 
0.95 
0.96 
0.95 
0.95 


0.96 
0.96 
0.97 
0.98 
0.93 



Table 5: Inclusive dijet cross-sections da j dif . Other details as in the caption to 
TableUl 
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lo Sio (0 bin 


da/d\og 10 (f) 
(Pb) 


"stat 


$syst 


#ES 


Cqed 


Chadr ' C*Z 


-2 ... 
-1.6 ... 


-1.6 
-1.45 


62.63 
143.3 


±0.91 

±2.1 


+0.81 
-0.86 

+2.5 
-2.5 


+3.27 
-2.83 

+7.1 
-6.6 


0.97 
0.96 


0.95 

0.95 


-1.45 ... 


-1.3 


143.0 


±2.1 


+1.0 
-0.8 


+7.3 
-6.1 


0.95 


0.96 


-1.3 ... 


-1.1 


109.9 


±1.5 


+2.9 
-3.0 


+4.9 
-4.9 


0.95 


0.96 


-1.1 ... 





17.40 


±0.24 


+0.11 
-0.13 


+0.65 
-0.63 


0.94 


0.98 



Table 6: Inclusive dijet cross-sections do j d log i0 (£). Other details as in the 
caption to TableUi 
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l°gio(0 bin 




da/d\og w (£) 














(pb) 


Ostat 


^syst 




Cqed 


Chadr " Cz° 


125 < Q 2 < 250 GeV 2 


-2.1 ... 


-1.65 


30.34 


±0.60 


+0.51 
-0.55 


+1.58 
-1.38 


0.97 


0.95 


-1.65 ... 


-1.5 


76.2 


±1.5 


+1.0 
-0.9 


+4.5 
-4.0 


0.96 


0.94 


-1.5 ... 


-1.3 


63.2 


±1.2 


+1.1 
-1.0 


+3.8 
-3.3 


0.97 


0.97 


-1.3 ... 


-0.4 


11.53 


±0.22 


+0.31 
-0.33 


+0.60 
-0.61 


0.97 


0.94 


250 < Q 2 < 500 GeV 2 


-2 ... 


-1.55 


19.93 


±0.52 


+0.29 
-0.27 


+0.82 
-0.71 


0.96 


0.95 


-1.55 ... 


-1.4 


48.2 


±1.3 


+2.1 
-2.1 


+2.3 
-2.2 


0.94 


0.96 


-1.4 ... 


-1.25 


42.8 


±1.2 


+0.5 
-0.4 


+2.1 
-1.9 


0.95 


0.96 


-1.25 ... 


-0.4 


8.56 


±0.21 


+0.13 
-0.14 


+0.40 
-0.41 


0.95 


0.95 


500 < Q 2 < 1000 GeV 2 


-1.9 ... 


-1.45 


8.21 


±0.29 


+0.22 
-0.19 


+0.29 
-0.27 


0.94 


0.94 


-1.45 ... 


-1.3 


28.32 


±0.93 


+0.74 
-0.72 


+ 1.10 
-0.76 


0.94 


0.96 


-1.3 ... 


-1.15 


29.11 


±0.93 


+1.23 
-1.27 


+1.22 
-1.02 


0.94 


0.97 


-1.15 ... 


-0.4 


6.00 


±0.17 


+0.07 
-0.08 


+0.26 
-0.25 


0.95 


0.97 


1000 < Q 2 < 2000 GeV 2 


-1.7 ... 


-1.25 


4.91 


±0.22 


+0.16 
-0.16 


+0.13 
-0.15 


0.94 


0.95 


-1.25 ... 


-1.15 


15.52 


±0.80 


+0.64 
-0.63 


+0.80 
-0.27 


0.93 


0.98 


-1.15 ... 


-1 


16.87 


±0.68 


+0.10 
-0.15 


+0.41 
-0.58 


0.94 


0.98 


-1 ... 


-0.25 


2.97 


±0.12 


+0.06 
-0.06 


+0.11 
-0.10 


0.93 


1.00 


2000 < Q 2 < 5000 GeV 2 


-1.5 ... 


-1 


3.11 


±0.16 


+0.07 
-0.08 


+0.04 
-0.06 


0.93 


1.03 


-1 ... 


-0.85 


9.07 


±0.48 


+0.26 
-0.28 


+0.12 
-0.18 


0.92 


1.03 


-0.85 ... 


-0.2 


2.54 


±0.12 


+0.09 
-0.09 


+0.08 
-0.07 


0.93 


1.04 


5000 < Q 2 < 20000 GeV 2 


-1.1 ... 


-0.75 


0.865 


±0.099 


+0.065 
-0.065 


+0.011 
-0.013 


0.94 


1.10 


-0.75 . . . 


-0.55 


2.85 


±0.23 


+0.10 
-0.10 


+0.02 
-0.03 


0.91 


1.08 


-0.55 ... 





0.794 


±0.071 


+0.092 
-0.090 


+0.023 
-0.020 


0.92 


1.10 



Table 7: Inclusive dijet cross-sections da / d\og 1Q (£) in different regions of Q 2 . 
Other details as in the caption to Tabled 
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pJ et Kir, 
Ulrp g um 




d(T j Ulrp g 












{\je\ ) 




(^pD/ Vjre V ) 


Ostat 


^syst 


Ses 


L'QED 










125 < Q 2 


< 250 GeV 2 






8 ... 


15 


5.050 


±0.057 


+0.071 
-0.070 


+0.311 
-0.274 


0.97 


0.95 


15 ... 


22 


1.385 


±0.028 


+0.037 
-0.038 


+0.063 
-0.062 


0.97 


0.96 


22 ... 


30 


0.292 


±0.012 


+0.012 
-0.013 


+0.009 
-0.010 


0.97 


0.96 


30 ... 


60 


0.0241 


±0.0016 


+0.0009 
-0.0008 


+0.0011 
-0.0010 


0.97 


0.95 








250 < Q 2 


< 500 GeV 2 






8 ... 


15 


2.937 


±0.046 


+0.076 
-0.076 


+0.146 
-0.141 


0.95 


0.95 


15 ... 


22 


0.998 


±0.026 


+0.011 
-0.011 


+0.040 
-0.035 


0.95 


0.98 


22 ... 


30 


0.215 


±0.011 


+0.008 
-0.008 


+0.006 
-0.008 


0.96 


0.97 


30 ... 


60 


0.0195 


±0.0016 


+0.0015 
-0.0015 


+0.0010 
-0.0006 


0.93 


0.95 








500 < Q 2 


< 1000 GeV 2 






8 ... 


15 


1.502 


±0.031 


+0.055 
-0.054 


+0.064 
-0.052 


0.94 


0.95 


15 ... 


22 


0.629 


±0.019 


+0.008 
-0.009 


+0.023 
-0.021 


0.95 


0.99 


22 ... 


30 


0.1665 


±0.0089 


+0.0041 
-0.0041 


+0.0054 
-0.0040 


0.96 


0.98 


30 ... 


60 


0.0194 


±0.0015 


+0.0010 
-0.0010 


+0.0007 
-0.0010 


0.95 


0.99 








1000 < Q 2 


< 2000 GeV 2 






8 ... 


15 


0.701 


±0.020 


+0.017 
-0.017 


+0.025 
-0.022 


0.93 


0.95 


15 ... 


22 


0.352 


±0.014 


+0.012 
-0.013 


+0.011 
-0.009 


0.94 


1.01 


22 ... 


30 


0.0943 


±0.0064 


+0.0063 
-0.0063 


+0.0025 
-0.0026 


0.94 


1.02 


30 ... 


60 


0.0136 


±0.0012 


+0.0003 
-0.0004 


+0.0006 
-0.0007 


0.94 


1.04 








^ AAA , /\V. 

2000 < Q z 


< 5000 GeV 2 






8 ... 


16 


0.350 


±0.013 


+0.009 
-0.009 


+0.004 
-0.007 


0.92 


1.00 


16 ... 


28 


0.1191 


±0.0058 


+0.0023 
-0.0022 


+0.0030 
-0.0023 


0.93 


1.07 


28 ... 


60 


0.01040 


±0.00097 


+0.00053 
-0.00049 


+0.00044 
-0.00046 


0.94 


1.08 








5000 < Q 2 


< 20000 GeV 2 






8 ... 


16 


0.0995 


±0.0076 


+0.0092 
-0.0092 


+0.0012 
-0.0005 


0.93 


1.05 


16 ... 


28 


0.0354 


±0.0031 


+0.0023 
-0.0021 


+0.0003 
-0.0008 


0.89 


1.14 


28 ... 


60 


0.00368 


±0.00053 


+0.00015 
-0.00023 


+0.00016 
-0.00012 


0.95 


1.20 



Table 8 

details as in the caption to Tabled 
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Figure 1: Uncorrected data distributions for inclusive dijet production (dots). 
For comparison, the predictions of the ARIADNE (dashed histograms) and LEPTO 
(solid histograms) MC models are also included. 
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Figure 2: Comparison of uncorrected data (dots) and MC model predictions for 
distributions of \og 10 (£) in different regions of Q 2 . For comparison, the predictions 
of the ARIADNE (dashed histograms) and LEPTO (solid histograms) MC models 
are also included. 
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Figure 3: Comparison of uncorrected data (dots) and MC model predictions for 

distributions of E^ B in different regions of Q 2 . For comparison, the predictions 
of the ARIADNE (dashed histograms) and LEPTO (solid histograms) MC models 
are also included. 
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Figure 4: The fraction of gluon-induced events as a function of \og 10 (£) as 
predicted by the CTEQ6.6 PDFs in different regions of Q 2 . 
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Figure 5: The fraction of gluon-induced events as a function of E } ^ B as predicted 
by the CTEQ6.6 PDFs in different regions of Q 2 . 
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Figure 6: The relative CTEQ6.6 PDF uncertainty, the relative uncertainty due to 
missing higher orders estimated by a variation of fiR and the theoretical predictions 
from different PDF sets relative to those obtained with CTEQ6.6 as functions of 
log 10 (£) in different regions of Q 2 . 
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Figure 7: The statistical uncertainty, 5 s tat ; the correlated uncertainty associated 
with the energy scale of the jets, 5es, an d the quadratic sum of the correlated and 

uncorrelated, 5 syst , systematic uncertainties, J o| s + 6® yst , as functions of Q 2 . 
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Figure 8: The measured differential cross-sections (a) da/dQ 2 and (b) 
da/dxsj for inclusive dijet production with E^ B > 8 GeV, M# > 20 GeV and 

— 1< //lab < 2.5 (dots), in the kinematic range given by 0.2 < y < 0.6 and 
125 < Q 2 < 20000 GeV 2 . The inner error bars represent the statistical uncer- 
tainty. The outer error bars show the statistical and systematic uncertainties, not 
associated with the uncertainty on the absolute energy scale of the jets, added in 
quadrature. The shaded bands display the uncertainties due to the absolute energy 

— ~ 2 

scale of the jets. The NLO QCD calculations with n 2 R = Q 2 + £y B (solid lines), 

— ~ 2 

jj, 2 R = Q 2 (dashed lines) and n 2 R = E } ^ B (dotted lines), corrected for hadronisation 
effects and Z° exchange and using the CTEQ6.6 parameterisations of the proton 
PDFs, are also shown. The lower parts of the figures show the relative differences 

— ~ 2 

with respect to the NLO QCD calculations with \i 2 R = Q 2 + E^ B ■ The hatched 
bands display the total theoretical uncertainty. 



27 



ZEUS 



*t 10 
O 

is 

a, 

T5 



10 1 

o 

^0.2 

* o 
Sg 

■3 -0.2 



• ZEUS (374 pb" 1 ) 



NLO® C. . ® C„. 

hadr Z 




~ \/ / \ NLO uncertainty 

I i i i i I i i i i I i i i i 



10 



20 



30 



40 



50 



60 

E^(GeV) 




Figure 9: The measured differential cross-sections (a) da/ dE ] ^ B and (b) da/ dM^ 
for inclusive dijet production. Other details as in the caption to Fig. 
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Figure 10: T/ie measured differential cross-sections (a) da j drf and (b) 
da/d\og 10 (£) /or inclusive dijet production. In the last r]* bins the NLO QCD 
predictions are not plotted for reasons explained in the text. Other details as in the 
caption to Fig. 
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Figure 11: The measured differential cross-section do j d \o<g 10 (£) for inclusive 
dijet production in different regions of Q 2 . Other details as in the caption to Fig. 
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Figure 12: Relative differences between the measured differential cross- 
sections da I dlog 10 (£) presented in Fig. [71] and the NLO QCD calculations with 
— ~ 2 

I^r = Q 2 + ^tb ■ Other details as in the caption to Fig. [3 
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Figure 13: The measured differential cross-section da/dE 3 ^ B for inclusive dijet 
production in different regions of Q 2 . Other details as in the caption to Fig. 
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Figure 14: Relative differences between the measured differential cross- 
sections da/dEj,-Q presented in Fig. [73] and the NLO QCD calculations with 
— ~ 2 ' 

t l R = Q 2 + Etb ■ Other details as in the caption to Fig. 
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